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We study a one-dimensional multicomponent anyon model that reduces to a multicomponent 
Lieb-Liniger gas of impenetrable bosons (Tonks-Girardeau gas) for vanishing statistics parameter. 
At fixed component densities, the coordinate Bethe ansatz gives a family of quantum phase transi- 
tions at special values of the statistics parameter. We show that the ground state energy changes 
extensively between different phases. Special regimes are studied and a general classification for the 
transition points is given. An interpretation in terms of statistics of composite particles is proposed. 
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Introduction — Anyons play an important role in topo- 
logical quantum computational They are the natu- 
ral quasiparticle excitations of fractional quantum Hall 
statesP While the fractional quantum Hall effect occurs 
in a two-dimensional electron liquid, its e dge excitations 
are described by one-dimensional anyonsPEI These edge 
states are essential for performing quant um c omputation 
operations by topological quantum gatesPZl 

A model for one-dimensional anyons is the Lieb-Liniger 
modeP with fractional exchange statistics.^ An impor- 
tant feature of this model is that it is solvable by the 
Bethe ansatzP^I I n this model the statistics parame- 
ter (exchange phase) does not modify the bulk ground 
state energy in large systems (correction terms vanish 
in the thermodynamic limit): 11 However, the effect of 
anyon statistics on other nonlocal physical quantities 
are more prominent. For example, momentum distri- 
butions, reduced density matrices, and other correlation 
functions have a nontrivial dependence on the stati stics 
parameterl^^H Experiments involving photoniP^U anc j 
ultracold atoms^ that have anyonic exchange statistics 
in one-dimension are accessible with current technology. 

In this paper, we generalize the anyonic Lieb-Liniger 
model to the multicomponent case. As known for the 
integrable multicomponent Bose gaspSHUl the additional 
internal degrees of freedom lead to diverse physical 
effect s paESl f n particular, we consider a gas of two mutu- 
ally impenetrable anyonic species. This exactly solvable 
model may prove useful in studies of the edge states of bi- 
layer fractional quantum Hall systems PaSa We encounter 
an interesting quantum phase transition^ (QPT) that 
is controlled by the statistics parameters of the model. 
This QPT distinguishes between two phases with ground 
states that are described by either one or two Dirac seas 
(one-dimensional Fermi spheres). The ground state en- 
ergies of these phases can therefore differ significantly by 
an extensive amount. 

We first motivate this work by describing a similar 
QPT that is controlled by the dynamical interparticle 
coupling in a multicomponent boson gas. We then pro- 
ceed with a formal description of our multicomponent 
anyonic model. This model is formulated in terms of 
quantum fields satisfying generalized commutation rela- 
tions. Next, we calculate exact eigenfunctions and the 



energy spectrum of this theory. We finally show how the 
statistics parameters determine whether the ground state 
is described by one or more Dirac seas. This result allows 
us to construct a phase diagram for the ground state in 
the space of all possible statistics parameters. 

Bosons — Let us consider a generalized two-component 
boson gas in the impenetrable or Tonks-Girardeau (TG) 
limitP^Hereafter, we sometimes use a pseudospin picture 
and refer to these two components as different spin-1/2 
projections. The relevant quantum Hamiltonian is a mul- 
ticomponent version of the Lieb-Liniger model 

2 ,L 



Hb = E / [d x 4,i{x))[dMx)} 

a=l J ° 



+ ^2c a p<j)i{x)(j)' !< l3 {x)(f) l3 (x)<f> a (x)dx. (1) 
0=1 

The fields 4> a (x) are independent canonical Bose fields 
satisfying periodic boundary conditions in a box of length 
L. The field index a € {1,2} labels different par- 
ticle species. Energy eigenstates of this Hamiltonian 
are simultaneous eigenstates of the number operator 
Na = Jq 4>a( x )4 > a(x) dx. Thus, the number of type-ev 
particles N a are good quantum numbers. 

We look at two extreme limits of the interspecies cou- 
pling. One is the decoupled limit in which the matrix 
c = ( cW C22 ) i s diagonal. That is, C12 = c 2 i = while 
both en and c 2 2 tend to infinity. This limit describes 
two non-interacting TG gases. The quantum mechani- 
cal eigenfunctions of the model are Bethe wavefunctions. 
The impenetrable nature of the gas implies that this 
wavefunction vanishes when the coordinates of particles 
of the same spin coincide. Thus, the wavefunction xb is 
a product of two symmetrized Slater determinants: 



det e iAl 'J Zlfc 



X^(zi,z 2 |Ai, A 2 ) = J == 



|dete iA2 ' mZ2 H , , 
xJ t== (2) 



Each N a x N a determinant depends only on the coordi- 
nates z Q = {z a ^} of type-a particles, i € {1,...,A Q }. 
They are characterized by independent sets of N a spec- 
tral parameters or momenta X a = {A Qi i}. The momenta 
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within each set are distinct, otherwise the wavefunc- 
tion ([2| vanishes (Pauli exclusion).^ Hence, the ground 
state is described by two Dirac seas with Fermi lev- 
els q a — (N a — l)n/L. In the thermodynamic limit 
where particle numbers N a tend to infinity at fixed lin- 
ear density D a , the ground state energy per unit length 
is Etf/L = (D\ + D^)tt 2 /3. 

The opposite limit is of mutual repulsion where all el- 
ements of c tend to infinity.^ In this example both com- 
ponents are coupled dynamically. The eigenfunctions in 
this case are 



Xb(z,<t n |A) = A(cr N ) 



|dete iA ^ fe | 
7WT~ 



(3) 



In this formula, the vectors z and <tn give the position Zi 
and spin <Ji of the i th boson. The vector <tn has N a com- 
ponents equal to a. The factor A(<x N ) gives the overall 
symmetry of the wavefunction under particle exchanges. 
It is equal to +1 when coordinates of particles with iden- 
tical spins are interchanged. Otherwise, when coordi- 
nates of different spins are exchanged it is equal to +1 
or —1. A choice of ±1 will correspond to a compatible 
irreducible repre senta tion (Young tableau) of the sym- 
metric group 6*2 ! 46 l 47 l Meanwhile, the wavefunction ^ 
is characterized by a single set A = {Xj} of N distinct 
momenta. The ground state of this two-component gas 
is thus described by a single Dirac sea with Fermi level 
q = (N— l)7r/£P^In the thermodynamic limit with total 
particle density D = D% + D2, the ground state energy 
per unit length is E\/L = D 3 tt 2 /3. This value is equal 
to the ground state energy density of a single component 
TG gas and fully-polarized free Fermi gas with the same 
total particle number? 32 ! 49 ! 

The ground states of these boson gases are fundamen- 
tally different. One is characterized by a single Dirac sea 
and the other by two Dirac seas. The ground state en- 
ergy of the mutually repulsive gas is much higher than the 
non-interacting example because [D\ + D2) 3 > Df + D\ 
(equal only in the "ferromagnetic" case when one of D a — 
0). Restructuring of the ground state between these two 
limits is controlled dynamically, that is, through the in- 
terspecies coupling coefficients. In the following section 
we seek to find a similar transition that is generated sta- 
tistically by generalized exchange phases between parti- 
cles. 

Anyons — Our one-dimensional anyonic model consists 
of two anyon fields ^(x) with a G {1,2}. These fields 
satisfy the exchange relations (no sums involved) 

^ a {x)^p{y) = U a[j (x - y)V(i(v)* a {x), 

*U*)*)>(v) - U afj (x - y)*l(y)¥ a (x), (4) 

¥a(aO*£(iO = Kp( x - tf)*J(»)*a(a:) + S aP 6(x - y), 

where U afj (x) = e 2 ™<*e< x \ Here K = Zl ) is a real 
symmetric matrix and e(x) is the sign of x with e(0) = 0. 
All cases of interest are contained in the parameter space 
K Q( g € [0, 1]. The exchange relations Q therefore depend 



on continuous parameters that interpolate between com- 
muting (bosonic) and anticommuting (fermionic) fields 
at different points in space and equal times. At the same 
space-time point we can choose to recover either bosonic 
fields or fermionic fields.^ The commutation relations 
here reduce to bosonic ones at the same point. 

The exchange phase of a localized bunch of N anyons 
can be computed from these exchange relations. Let 
&(Z) be a composite operator &(Z) = .(zi) 
for some configuration of anyons with positions z and 
spins er N . Further, let all {z,} be in the neighborhood of 
Z and all {yi} be in the neighborhood of Y so that all 
Zi > yj for some Z > Y. Hence, the exchange relation 
satisfied by a composite particle resulting from the fusion 
of N anyons is, for Z > Y, 

$ t (Z)$ t (F) = e 27ri (KiiW 1 2 +2 Kl2 JViA' 2 +K 22 Ar 2 2 )^ ( t(y)4 ) t^) 

^' (5) 
Now, the quantum Hamiltonian we consider is an any- 
onic version of the mutually repulsive boson model ([I]) 
with all c a p —¥ 00. It is given by 



^ = E [ L [d x *t(x))[d x * a (x)] 
a =i Jo 



+ cJ2H(x)*Ux)*p(x)y<x(x)dx. (6) 

0=1 

where c is a coupling constant that we will send to +00. 
In this section, we solve for the eigenfunctions and spec- 
trum of this model by the coordinate Bcthe ansatz. 

The component number density operators p a (x) = 
^ r J e (x)^ r a (x) obey the usual commutation relations 

[* a (a0, Pp(v)\ = S aP S(x - y)^fi(v), 
{^ a (x),pp(y)} = -6 aP 5(x-y)*l(y). (7) 

We find that the numbers N a of type-a anyons are indi- 
vidually conserved. 

We stress that this anyonic model (pi is related, but 
not equivalent, to the bosonic model (flT! To see this, we 
use a generalized Jordan- Wigner transformation 

^( I ) = e 2,ri ^ K °^»"^ Wdz * a (i), 

4(2/) = *tfo)e-™£^/oW«)«* ( 8 ) 

to obtain independent canonical Bose fields 4> a (x). Ap- 
plying the inverse transformation to the anyonic model 
([6| results in a gauge-transformed Lieb-Liniger model 

2 „L 

K=/Z [D a {x)^ a {x)}\v a {x)^ a {x)} 

a = l J ° 

2 

+ c ^4'ai. x )4'p(x)</>p(x)(l) a (x)dx, (9) 
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with covariant derivative T> a (x) = d x — 2ni ^ 7 K Q7 p 7 (x). 

Following the treatment of the single component any- 
onic caseP we map the Hamiltonian 01 to a many-body 
quantum mechanical problem. We let the vectors z and 
<tn contain the positions Zi and spin Gi of the i th anyon, 
as done with the mutually repulsive boson gas pi). We 
introduce the number eigenstates (TV = Ni + N 2 ) 



\N U N 2 ) 



( CT N> 



N 

j d jv z X (z, < r N )n^(^)io), (io) 



where the anyonic vacuum state satisfies ^ a (zi)\0) = 
0. The sum over (o"n) denotes a sum over all con- 
figurations with N a type-a anyons. This construction 
imposes a certain anyonic symmetry on the wavefunc- 
tion x- That is, if Vij is a particle exchange op- 
erator such that Vijxi- ■ ■ z % ■ ■ ■ z j • ••■>••• ■ ■ ■ &j ■ ■ ■ ) = 
X(. ..Zj... Zi (Xj ... <Ji ... ), then 

^ — e 27ri E? =i+ i «<r i n e fe _2; ') e - 27ri ECi+i K « j <'i e ( z 3- Z i)'p. 

'fa 

The associated many-body Hamiltonian acting on the 
wavefunction %(z, <tn) is 



/V 



2c 



E 



- 2T,-). 



(12) 



In the Tonks-Girardeau limit c — > oo the wavefunction x 
vanishes at the coinciding points Zi — Zj (i =/= j). With 



the symmetry relations (111, we are able to write the any- 



onic wavefunction % as a product of a multicomponcnt 
boson wavefunction xb(z, <tn) and a position dependent 
anyonic phase: 



x(z,<r N ) = XB(z,crN)e 



(13) 



The symmetry of the anyonic wavefunction x between 
an interchange of anyons of different spins is not fixed 
by the anyonic commutation relations Q. Like in the 
bosonic case ([3]), the wavefunction xb m Eq. (13) con- 
tains a factor corresponding to the appropriate irre- 
ducible representation of the symmetric group S 2 . The 
choice of Young tableau is related to the particular set 
of fusion rules that we can assign to the anyons of the 
model. For simplicity we consider the completely sym- 
metric Young tableau of N = N\ + N 2 boxes for the 
boson wavefunction. That is, we assume that xb in 



Eq. ( 13 ) is also symmetric between interchanges of dif- 



ferent types of particles. The antisymmetric case may be 
obtained from the symmetric case by the transformation 

«12 -> «12 + 1/2. 

With these considerations, the anyonic wavefunction 
can be explicitly expressed as 



X(z,ctn|A) 



dot c 1 



cm 



A] , 

A 2 -* 



-Qch- 



XL I 2k 

FIG. 1. (Color online) Ground state occupation of momen- 
tum vacancies for two type-1 and four type-2 impenetrable 
anyons. Filled vacancies are shaded and the maximum oc- 
cupancy is one. Squares and circles refer to different sets of 
vacancies. Each set is shifted from its bosonic limiting value 
by a statistics-dependent term A a . There are two Dirac seas 
when the shifts are unequal Ai 7^ A2 (top) and one Dirac sea 
when they are equal Ai = A2 (bottom). 



where the set of N momenta {A m } are determined by 
the boundary conditions imposed on x- We will find it 
convenient to relabel the momenta {A m } — ¥ {X a ,i} with 
Greek index a G {1, 2} and Latin index i £ {1, . . . , iV Q }. 
Hence, a state with this wavefunction has total momen- 
tum P = E Q Ei V* and energy E = J2 a E* 

The boundary conditions on the anyonic wavefunction 
may be chosen to be periodic or twisted. If we choose 
the twisted boundary condition to cancel the anyonic 
phase e~ m ^ ,i<j K °i"j e ( Zi ~ z 3) 5 W e recover the energy spec- 
trum of the mutually repulsive boson model. This phe- 
nomenon is also observed in the single component caseP^ 
That is, certain twisted boundary conditions on the any- 
onic wavefunction results in a bosonic wavefunction and 
vice versa. This correspondence is related to the Jordan- 
Wigner transformation Q connecting the anyonic model 
Q and gauge-transformed bosonic model ffim. 

We now focus on the situation where the anyonic wave- 
function x satisfies periodic bound ary co nditions. As dis- 
cussed for the one component caseP*E3 periodic bound- 
ary conditions on x enter in a nontrivial way. Suppressing 
component labels, the periodic boundary conditions we 
impose on x( z ii z 2, ■ ■ ■ , z n) ar e 



X(0,z 2 ,z 3 ,...) =x(L,z 2 ,z 3 ,...) 
X (zi, 0, z 3 ,...) = e - 4 ™*^ X {zi,L, z 3 , ... ) 
X (zi,z 2 ,0, . . . ) - e-M^an+S.)^,^^, . . . ) 
X (z k = 0) = e -**X£l «*>- ix ( Zk = L). (15) 

In logarithmic form, the Bcthe equations for the anyonic 
momenta are therefore 



■ + Tr(N - 1) + 5 a 



(16) 



where {n a i ] is a set of integers. Here ir(N — 1) is the ac- 
cumulated dynamical scattering phase due to the anyon 
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hard cores. The phase shift S a is 



S a = 27T| ■ /,-,,,, + 2 «. 

7=1 



(17) 



This quantity represents the total phase shift due to the 
exchange statistics. We can also write the Bethe solutions 
(fl6| as 



A, 



2tt 



(m ai + A Q ), 



(18) 



where {m Qj i} is a set of half-odd integers for even N (or 
integers for odd N). The shift A a is defined as 



fraC (fe)' 



fracf 
I \2ttJ 



if cfjj) € [0,1/2], 



1. if frac 



(|) ed/2,1). 



(18) 



where frac(x) is x minus the integer part of x. With this 
definition, A Q is restricted to the interval (—1/2, 1/2]. 
Quantum phase transition — The momenta that label 



the anyonic wavefunction ( 14 ) are different from each 
other. Otherwise the Slater determinant is zero and the 
wavefunction vanishes. Thus, when the shifts Jl9|) are 



equal Ai = A2, the ground state is described by a single 
filled Dirac sea. That is, N momentum vacancies {A m } 
with the smallest squares are occupied (see Fig. [T]) . On 
the other hand, when Ai 7^ A 2 the set of momentum va- 
cancies split into two: {A^} and {\2.j}- Two Dirac seas 
are therefore formed in the ground state. Each is filled 
with N a quasiparticles to minimize the ground state en- 
ergy. Since the shifts A a depend on the parameters K a p , 
we may think of the transition between the single and 
double Dirac sea ground states as QPTs controlled by 
K a p. Hereafter, we will refer to the single Dirac sea phase 
as degenerate and double sea phase as nondegenerate. 

The ground state of the model generally has nonzero 
momentum as the anyonic statistics break spatial par- 
ity symmetry. In both phases, the ground state acquires 
a total momentum P = 2ir(D 1 A 1 + D 2 A 2 ) for N X ,N 2 
even. Thus, the statistics parameters modify the ground 
state energy by a constant (non-extensive) amount. How- 
ever, the ground state energy of these phases are signifi- 
cantly different from each other. The exact ground state 
energy density for the degenerate phase is 



L 



TV 

y 



(£>! + D 2 ) 



(12A 2 -!)(£!+ ,D 2 ) 
L 2 



, (20) 



where A = Ai = A 2 . Meanwhile, for the nondegenerate 
phase it is equal to 



El 
L 



D\ + Dl 



(12A 2 - 1)1?! + (12A 2 - 1)D 2 
L 2 



for Ni , N 2 even. For an anyonic gas with fixed component 
densities, the degenerate phase has a higher ground state 




FIG. 2. (Color online) Phase diagram for a gas of two type-1 
and four type-2 impenetrable anyons. The (parallel) planes 
represent values of K a p for which the ground state is charac- 
terized by a single Dirac sea. 



energy than the nondegenerate phase (equal only in the 
ferromagnetic situation). This difference persists even 
in the thermodynamic limit of large systems and fixed 
densities. 

Let us further examine the features of this quantum 
phase transition by considering some special cases of in- 
terest. The first is the limit of a single-component anyon 
gas. That is, all statistical parameters are equal to each 
other K a p = k. In this case the degeneracy condition 
Ai = A 2 is always satisfied and the ground state is de- 
scribed by a single Dirac sea. This case includes the mu- 
tually repulsive boson model (pi) (k = 0) and polarized 
pseudofermionic model (n = 1/2). 

Second, we take the matrix k to be diagonal. Here, 
the accumulated statistical phase 5 a only depends on the 
number of a-type anyons. Thus, the only coupling be- 
tween components is dynamic in nature and due to hard 
core repulsion. The condition Ai = A 2 for a degenerate 
phase reduces to 



Kll(N! - 1) - K 2 2(N 2 - I) =p, 



(22) 



where p is an integer. When the components have iden- 
tical anyonic statistics K\\ — k 22 the ground state is 
characterized by a single Dirac sea in the balanced case 
iVi = N 2 . 

Third, we consider the case in which the only nonzero 
elements of k are off-diagonal. We have two mutually 
impenetrable boson gases, but with fractional exchange 
statistics between particles of different types. The ground 
state of the system becomes degenerate when 



Kia(JVa-JVi) =p. 



(23) 



This condition is fulfilled automatically in the case of a 
balanced mixture N± — N 2 . Otherwise, degeneracy will 
generally not happen for irrational k± 2 - 

Next, we consider a mixture of particles and anti- 
particles. In this case the statistics parameters satisfy 2 
Kn = k 22 — —K\2. Substituting these conditions into the 
rule Q shows that an equal mixture of particles and anti- 
particles (Ni = N 2 ) has no overall exchange phase. This 
is the expected result because the vacuum obtained from 
pair annihilation has no phase. For an arbitrary mixture 
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of particles and anti-particles the degeneracy condition 
becomes 



2K U {Nt - N 2 ) = p. 



(24) 



Hence, the degenerate ground state always exists for the 
balanced case Ni = N 2 and also for the pseudofermionic 
situation k±i — 1/2. 

In the most general case, the degeneracy condition 
Ai = A2 is satisfied when 

(JVi - 1)ku + (N 2 - JVi)«ia - (N 2 - 1)k 22 = p. (25) 

The single Dirac sea phase occurs on a family of planes 
in n a p space. Let us write the vectors K = k 12 , K22) 
and n = (N\ — 1, N2 — Ni, —N 2 + 1) in the same orthonor- 
mal basis. The condition (25) then becomes 



K-p 



n 



•n = 0, 



(26) 



which is a set of standard equations of planes labeled by 
the integer p. The distance d between neighboring planes 
is d = [(iVi - l) 2 + (N 2 - m) 2 + (N 2 - l) 2 ]- 1 / 2 . Thus, 
in the cube K a p € [0, 1], the generic situation for a finite 
number of particles is the nondegenerate double Dirac 
sea phase. A sample phase diagram is given in Fig. [2j 

The set of planes given by Eq. ( 25 ) denote the values 
of the statistical parameters K Q( g and component num- 
bers N a for which the ground state is characterized by 
a single Dirac sea (degenerate phase). For fixed compo- 
nent densities, a parametric change of the statistics pa- 
rameters K, a p{t) will cross these quantum critical planes 
at certain points (the parameter t describes a curve in 



n a j3 space). These quantum critical points define phases 
where the ground state energy is much larger than the 
sum of ground state energies of each component. 

Concluding remarks — We discovered that anyonic 
statistics in a multicomponent Tonks- Girardeau gas can 
lead to the formation of a new ground state with in- 
teresting physical features. The Bethe momenta in this 
state form two independent sets of vacancies. Each set is 
therefore filled independently with quasiparticles to form 
independent Dirac seas in the ground state. This state 
is remarkably different from the ground state of a mix- 
ture of impenetrable bosons, which has only one Dirac 
sea. For example, the formation of multiple Dirac seas 
significantly reduces the energy of the system and favors 
a balanced mixture of components (zero magnetization) . 
This effect competes with external chemical potentials 
that tend to drive the system toward homogeneity (fer- 
romagnetism) . 

Moreover, the quantum phase transition described here 
is not observed in the single-component anyonic Lieb- 
Liniger model. Unlike in the single-component case, the 
bulk ground state energy in the multicomponent model 
depends significantly on the statistics parameters even 
for large systems. 

The formalism presented here extends readily to any- 
onic TG gases with more than two components. In 



particular, the statistics-dependent phase shift 5 a (17) 
generalizes to S a = —2nn aa + 2ttJ2 7= i K a7 iV 7 for an M- 
component gas of impenetrable anyons. 
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